TDAE-C6Q: Can a Mott insulator be a ferromagnet? 
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Motivated by the structure of TDAE-Ceo, we derive a multicomponent superexchange Hamiltonian 
for the spins and orbital ("isospin") degrees of freedom in the Mott-insulating phase. We explore 
its phase diagram and identify points of special interest: an SU(4) antiferromagnet (solved by 
Sutherland) and two ferromagnet®antiferromagnet points where the ground state is known. For the 
ferromagnetic regime, we apply interchain mean field theory and derive an expression for the Curie 
temperature where spin ordering occurs and a lower Neel temperature for a conjectured isospin 
ordering. 
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The fullerene compound TDAE-Ceo, where Cqq is buck- 
minstcrfullerene and TDAE is tetrakis(dimethylamino)- 
ethylene C2N4(CH3)g, exhibits ferromagnetism at T c ea 
16°K QJ^]. The striking aspects of this discovery are (i) 
the magnitude of T c - relatively large for a material with 
no transition metals - and (ii) its nonmetallic conductiv- 
ity, suggestive of Mott-Hubbard localization Q . 

ESR studies H show that TDAE donates an electron 
to Ceo- Furthermore, no ESR signature of TDAE+ is 
observed, suggesting that the TDAE radical spins are 
somehow paired. The monoclinic structure makes for 
a relatively short inter-Cgo separation along the c-axis. 
We are then led to consider a model of Cg chains whose 
conduction electrons interact via superexchange. 

Superexchange in a one-band model is always antiferro- 
magnetic, since the intermediate state | Tl ) is a spin 
singlet. The molecular degeneracy of the t\ u Ceo LUMO 
leads to interesting possibilities not realized in a orbitally 
nondegenerate model. Indeed, Seshadri et al. || have 
discussed how ferromagnetism naturally arises via su- 
perexchange through intermediate states with a negative 
singlet-triplet splitting (e.g. Hund's rule) in Cqq. Moti- 
vated by Ref. || , we introduce here what we believe to be 
a 'minimal model', based on the structure of TDAE-Ceo, 
which leads to insulating ferromagnetic behavior. 

The main result of our Letter is the full multicompo- 
nent superexchange Hamiltonian for the Mott insulator, 
an analysis of its phase diagram and ground states, and 
a discussion of three-dimensional ordering in this quasi- 
one-dimensional system. The model is characterized by 
three coupling constants, which represent the superex- 
change interaction strength through the three different 
intermediate states of CgQ . In the limit where all the 
intermediate states are degenerate, we obtain an SU(4) 
Heisenberg antiferromagnet in its fundamental represen- 
tation, solved by Sutherland [|| using Bethe's Ansatz. 
In the case of large singlet-triplet splitting, we are able 
to prove using Marshall's theorem J?J that the ground 
state is fully polarized in the spin variables and is given 
by Bethe's wavefunction for the S = \ antiferromagnet 



in the isospin variables. This latter case, which we think 
relevant to TDAE-Ceo, is a rare example where ferromag- 
netism can be proven in a microscopic model of localized 
electrons ||. The three-dimensional susceptibilites can 
be expressed in terms of the exact one-dimensional func- 
tions following Scalapino et al. || . For the ferromagnetic 
model we find the spin Curie temperature and the orbital 
Neel temperature as a function of interchain couplings. 
Full details of the calculations below will be soon avail- 
able in a longer paper [jl8| . 

The Hopping Model: We consider tight binding hopping 
on a lattice of Ceo molecules with a filling of one electron 
per site. In general, a tetragonal or monoclinic crys- 
talline symmetry will resolve the triply degenerate t\ u 
orbital into three distinct levels. Better details could be 
obtained ab-initio once the precise structure and orienta- 
tions of the Ceo molecules arc experimentally ascertained. 
In our model, we shall retain only what we believe may 
be the essential microscopic physics undelying the ferro- 
magnetism in TDAE-Ceo^ (a) The hopping is quasi-one 
dimensional along the c-axis. (b) We assume that the 
crystal field resolves the t\ u orbital triplet into a lower 
doublet [l — ±) and a higher singlet I = at higher 
energy, as if the crystal fields are cylindrically symmet- 
ric about an axis which pierces the center a pentagonal 
face of Ceo- (c) Hopping along the chains is assumed to 
preserve the orbital magnetization I. 

Thus we investigate the Hamiltonian H = Hf wp + Ti.^ op + 
Ttion, where 

n L P = E (4(*te«r(« + c ) + H - c -) 

W hop = "I E *«'(*J-) (4r(0c,'„(» + S ±) + H - C -) 

1,1' , V 

W ion = E^ \ A d)) ( A (i)\ ■ (!) 
i,A 

Here cJ (i) creates, at site i, an electron of spin polar- 
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ization cr =],{ and "isospin" I = +, — . c and S± denote 
nearest neighbor lattice vectors in the c direction and 
the a-b plane respectively, tw are the hopping matrix 
elements between orbitals I and I' on neighboring chains. 

7iion is the interaction Hamiltonian which discourages 
multiple electron occupancy on any Cgo molecule. It is 
parametrized by pseudopotentials ua which correspond 
to the following multiplets, 

ui : c t +T c t „ T | ) , + ct +l c -T) I ) < c V-i I ) 

u 2 : c^ T c^ x | 0} ,cl T cl i | 0) (2) 

The relations between ua and the isotropic pseudopoten- 
tials ul of angular momenta L are: u\ = u\, u~2 = 112, 
but Mo = |mo + due to projecting out the "Z = 0" or- 
bital state. Thus, while in an isotropic environment there 
might be pair binding (up < 0) due to electron-electron 
|yj and electron- vibron |L2|] interactions, it does not pre- 
clude a repulsive uq > in the monoclinic crystal field 
environment. This may help to explain why TDAE-Cgo 
is not a CDW, nor a superconductor as is AsCqq. 

Multicomponent Superexchange Hamiltonian. Experi- 
ments have shown that TDAE-C60 is insulating at low 
temperatures, consistent with the existence of a gap to 
charge fluctuations (i.e. all u~a > 0) j|. The low-lying 
excitations can be described by a superexchange Hamil- 
tonian, formally obtained as a second order expansion 
in small t" ju. Since charge excitations are gapped, a 
renormalized version of the superexchange Hamiltonian 
is expected to describe the low energy excitations also for 

The zeroth order states of the superexchange Hamilto- 
nian are four singly occupied states enumerated by I, a. 
The operators which act on these states can be repre- 
sented by spin operators, S? = \ Ez )0 -, CT ' C L W T ao> c ia> (*) 

and "isospin" operators, I\ = \ YLi,y .0 c t(*) < T U' c v 
where r are the Pauli matrices. Taking into account the 
constraint ^2 a c\ t c a = 1, the 15 independent elements 
of the SU(4) generators Sp = cJ^Cg can be expressed in 
terms of the 15 operators {S^,r, S^I V }. 

One- Dimensional Limit. We first consider the purely 
one-dimensional case, where tn, = 0. There are three 
superexchange constants defined as 

J M = J^-L M = 0,1,2, (3) 
u M 

After a straigthforward calculation, we derive the general 
spin/isospin representation of the effective Hamiltonian 

as 

H = J2 i A S n ■ S n+1 + BI n -I n+1 +C (4) 

n 

+D S„ ■ S n +i I n ■ I n +1 + E S n ■ S n +i I^n+l + E) , 



where 

A = —\J\ + J2 + 2~Jo, B = I J\ — i Jo, 

C = J -J 2 , ^ = 2J!+2J , 

P = 4J 2 -4J , F = -§ J ± - \ J 2 - i J . (5) 

This model possesses a global SU(2) xU(l) symmetry, i.e. 
H commutes with J2 n S n and with J2 n ^n- Enlarged 
symmetries occur when J = J2, where the symmetry 
group is SU(2)xSU(2), and when J = J x = J 2 , where 
the symmetry group is SU(4). 

SU(4) Point - At the point u\ = U2 — uo = u (|J) ac- 
quires full SU(4) symmetry. For each c-chain the Hamil- 
tonian is 

^(4 ) = jEE^w s «( ,i + 1 )' ( 6 ) 

where J = 2{t^) 2 /u. 

The SU(P) Heisenberg antiferromagnet in the fundamen- 
tal representation has been solved by Sutherland for gen- 
eral P using Bethe's Ansatz ||. This model exhibits 
P — 1 gapless elementary excitation branches. We pre- 
sume, based on what happens in the SU(2) model p3[ , 
that for a chain of ./V sites where TV is an integer multi- 
ple of P, the ground state is an SU(P) singlet and the 
low-lying excitations transform according either to the 
singlet or the adjoint representation. This is essentially 
what happens in the fermion mean field (large P) theory 
of the SU(P) antiferromagnet Qfjl|[l6||. The mean field 
has four degenerate quarter-filled (fc F = jir) bands for 
P = 4. Although there it has no true long-ranged or- 
der, the spin and isospin susceptibilities diverges at the 
nesting wavevector 2fc F = ^7r, which describes a commen- 
surate spin density wave of period four. The period four 
arises because the spin chain is in its fundamental repre- 
sentation, and by '4-ality' one needs four sites to make a 
sine let 17|. The mean field theory also predicts a con- 
stant uniform (Pauli) susceptibility, and a linear specific 
heat as in a Fermi liquid |15| . 

Ferro-Antiferromagnetic points - Along the surface U2 = 
u , our Hamiltonian possesses an SU(2)xSU(2) symme- 
try. There are then two special limits in which we can 
determine the exact ground state, (i) The "F x A model" 
at Uq — > 00, with Jy = 2(t") 2 /ui, 

WfxA = — y^(SVi • S n+ i + §)(i — I n ■ I n +i) (7) 

n 

where the interactions are ferromagnetic in the spin chan- 
nel and antiferromagnetic in the isospin channel, and (ii) 
the "A x F model" for u\ — > 00, with J|| = 2(ill) 2 /u , 
and the roles of J and S interchanged. 
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It is possible to prove that the ground state of 7YfxA is 
the fully polarized ferromagnet | F ) s for the spin vari- 
ables, and Bethe's ground of the spin-half antiferromag- 
net for the isospin variables i.e. 



F ) s ® | Bethe } z 



(8) 



A corresponding result holds for Haxf, with spin and 
isospin variables exchanged. Due to the ££7(2) <g> SU{2) 
symmetry the total spin Stotj total isospin I tot , and their 
polarizations along the z axis (Ms and M/, respectively) 
are good quantum numbers. Following Lieb and Mattis' 
proof of the Marshall theorem for the Heisenberg model 
|l9| , we perform a it rotation about the z axis of the 
isospin operators on odd-numbered sites. The Hamilto- 
nian transforms into a non-positive ('negative semidefi- 
nite') operator in the product Ising basis 

n 

x(S n -S n+1 + l) = H' FxA (9) 



The accessibilty of all states within a given magnetiza- 
tion sector by repeated application of the Hamiltonian, 
implies (see Ref. [|l9|) that the ground state of H' FxA in 
the sector (Ms, Mj) = (0, 0) can be chosen to be positive 
definite in the sublattice-rotated Ising basis, i.e. it obeys 
Marshall's sign rule. Since the same Marshall signs hold 
for the state on the right hand side of Eq. B, which has 
Stot = \N, and I t ot = 0, the two sides of Eq. [s] have fi- 
nite overlap hence the same Stot and /tot ■ We are free to 
choose Ms = \N as a representative of the ground state 

manifold. Note that |tp^ xA ) is indeed an eigenstate of 
the spin triplet projection operator (S n ■ S n +i + |) with 
eigenvalue one. It follows from Eq. fj] that the isospin 
part of the wavefunction is the ground state of the spin- 
half antiferromagnetic Heisenberg chain, given by Bethe's 
Ansatz. 

Exact excitations of TCfxA within the isospin sector (re- 
taining full spin polarization) with dispersion ^ttJ\ sinfc| 
can be constructed as in Refs. |l3| , p0[| . The gapless ferro- 
magnetic magnons, which exist due to Goldstone's the- 
orem, can be approximated within the Single Mode Ap- 



proximation (SMA) : | k ) 
dispersion is 



s k I ^o FxA ) • The trial state 



uj(k) < 21n(2)J(l -cosfc) 



(10) 



from which we see that the ferromagnon bandwidth is 
decreased due to the antiferromagnetic nearest-neighbor 
isospin correlations, i.e. (j — J„ • I n +i) = ln(2). 

Classical Phase Diagram. The ground state depends on 
the dimcnsionlcss ratios uq/ui and U2/U1. The classical 
approximation (justified at S, I » 1) is given by mini- 
mizing the bond energies of Eq. (||) as function of vectors 



Si and Ij of magnitude \ . The results are plotted in Fig. 
%■ 

It is interesting to note that the SU(4) symmetry point 
is at the border of 4 distinct ordered phases of different 
symmetries, where the energy is degenerate along the 
lines (I n ■ I n +i) = —\ and (S„ ■ S n+ i) = \. The large 
degeneracy of the classical SU(4) model is reduced by 
quantum fluctuations. 

The classical regime of Heisenberg spin-ferromagnetism 
and isospin-antiferromagnetism extends throughout 
uo,U2 > U\, although quantum fluctuations break the 
SU(2) isospin symmetry away from the isotropic line 
uq = U2 (marked as a dashed line in Fig. |l|). 

3D Ordering in the F x A Model - As shown by 
Scalapino et al. one can treat the interchain inter- 
actions by mean field theory and thereby derive an ex- 
pression for the full susceptibility X a b(q±, Qz, in terms 
of (q z , , the susceptibility for the one-dimensional 
chains. The general result is 



X(q± , q z , uj) = [1 - J± (q± ) X 1D (q z , w)] 'x 1D (q z , w) 



(11) 



where J_l(<7_l) = J2s ± J±(^±) e ~ iq± S± IS trie spatial 
Fourier transform of the interchain coupling matrix. 
(Note that the quantities X, J±, and X 1D in Eq. |ll| are 
matrices.) This approximation also may be employed at 
finite temperature. 

Consider now the Fx A model discussed above. At fi- 
nite temperature T, long-ranged ferromagnetic order is 
destroyed and the global SU(2)xSU(2) symmetry is re- 
stored. The uniform susceptibility of the ferromagnetic 
chain is given by 



X'f(0,0;T) 



Jii 



24T 2 



as was first derived by Takahashi in Ref. |2l| 
Refs. 022|). 



(12) 



(see also 



For the antiferromagnetic susceptibility, we appeal to the 
bosonization results of Schulz and of Eggert and Affleck 
p3| , who have computed the dynamic susceptibility of 
the S = h antiferromagnetic Heisenberg chain. Perform- 
ing a Fourier transform of their result and taking the low 
frequency and wavevector limit near the antiferromag- 
netic point we obtain the staggered isospin susceptibility 



o_ 

irT 



(13) 



where ao ~ 4.44. 



For mixed interchain coupling operators e.g. O = S X I V , 
we may use the assumed independence of ferromagnetic 
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and antiferromagnetic magnons to obtain at low tem- 
peratures Xp A (7r, 0;T) ~ (J||T) -1 / 2 , which diverges even 
more slowly than Xa in the T — > limit. 

The interchain interaction is given by Jj_ = J\\((tj- ) 2 + 
(i^-) 2 )/4(tH) 2 , where t f - , t% are the transverse hopping 
integrals (see Ref. JOl). Thus, as the temperature is 
lowered, a transition from paramagnetic to ferromag- 
netic state should set in when J±X' F = 1. This yields 
Tc ~ y/ J\\ Jj_/24. The relation Tc oc ./Tjj Jj_ was also 
found by Scalapino et al. (ref. 0) in their studies of 
anisotropic Heisenberg magnets. It is conceivable that 
at still lower temperatures a Neel ordering of the isospin 
variables occurs at a Neel temperature Tn — 3a 2 , J±/ir. 

Experimental notes: (a) The lower isospin transition, 
to our knowledge, has not been resolved experimentally. 
Perhaps it is not very well separated from the ferromag- 
netic transition which would help explain the mysterious 
excessive entropy of transition found by Ref. (|] . 

(b) Alternatively, the isospin ordering might be pre- 
empted by a isospin-Peierls ordering (orbital dimeriza- 
tion) aided by the electron-phonon coupling. In that 
case, a signature for the isospin-Peierls effect should be 
present in X-ray scattering or in the phonon spectrum. 

(c) The role of a possible orientational disordered ground 
state j24| has not been considered here although it might 
help explain the observed weak ferromagnetism jlj]. In 
addition, Bloch's T 3 / 2 temperature dependence of the or- 
dered moment found in Ref. Q which holds upto T m Tc 
is hard to reconcile with quasi-one dimensionality where 
J± « T c . 

In summary, we have introduced a new model of quasi- 
one dimensional interacting electrons with doubly de- 
generate orbitals motivated by the structure of TDAE- 
Ceo- At occupancy of one electron per site, we obtain a 
Mott-insulator with muticomponent superexchange be- 
tween spins and isospins at neighboring sites. At spe- 
cial values of the interactions we identify exactly solvable 
points, including the SU(4) antiferromagnet, and spin- 
ferromagnet, isospin-antiferromagnet limit. The classical 
ground state diagram also contains a large region of spin 
ferromagnetism and orbital antiferromagnetism which we 
believe is relevant for TDAE-Cgo- A mean field analy- 
sis of the interchain coupling in this regime predicts two 
transition temperatures: ferromagnetic spin ordering at 
Tc oc •v/Jy Jj_, and orbital (isospin) antiferromagnetic or- 
dering at Tn oc Jj_. This lower transition, to our knowl- 
edge, has not yet been resolved experimentally. 
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S - AFM (Heisenberg) 
I - AFM (Ising) 



SU(2) x SU(2) 



-Q 



| FM ) s <g> | Bethe )_ 



S - FM (Heisenberg) 
I - AFM (Heisenberg) 




S - AFM (Heisenberg) 
I - FM (planar) 



S - AFM (Heisenberg) 
I - FM (Ising) 



Bcthe ) s <g> | FM ) 



FIG. 1. Ground state phase diagram for the Hamiltonian 
H of Eq. g. 
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